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General Structure of Relativistic Vector Condensation 

Francesco Sannino[] 

NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

We study relativistic massive vector condensation due to a non zero chemical potential associated 
to some of the global conserved charges of the theory. We show that the phase structure is very rich. 
More specifically there are three distinct phases depending on the value of one of the zero chemical 
potential vector self interaction terms. We also develop a formalism which enables us to investigate 
the vacuum structure and dispersion relations in the spontaneously broken phase of the theory. We 
show that in a certain limit of the couplings and for large chemical potential the theory is not stable. 
This limit, interestingly, corresponds to a gauge type limit often employed to economically describe 
the ordinary vector mesons self interactions in QCD. We finally indicate for which physical systems 
our analysis is relevant. 


I. INTRODUCTION 

Relativistic vector condensation has been proposed and studied in different realms of theoretical physics. However 
the condensation mechanism and the nature of the relativistic vector mesons themselves is quite different. Linde [Q, 
for example, proposed the condensation of the intermediate vector boson W in the presence of a superdense fermionic 
matter while Ambjprn and Olesen || investigated their condensation in presence of a high external magnetic field. 
Manton || and later on Hosotani [Q considered the extension of gauge theories in extra dimensions and suggested that 
when the extra dimensions are non simply connected the gauge fields might condense. Li in |^j has also explored a 
simple effective Lagrangian and the effects of vector condensation when the vectors live in extra space dimensions |p5| . 

Brown and Rlro || also suggested that the non gauge vectors fields such as the (quark) composite field p in QCD 
may become light and possibly condense in a high quark matter density and/or in hot QCD. Harada and Yamawaki’s 
jfj dynamical computations within the framework of the hidden local gauge Q symmetry support this picture. 

Rotational symmetry can also break in a color superconductor if two quarks of the same flavor gap. In this case 
the quarks must pair in a spin one state and a careful analysis has been performed in ||. Whether this gap occurs or 
not in practice is a dynamical issue recently investigated in @- 

We consider another type of condensation. If vectors themselves carry some global charges we can introduce a 
non zero chemical potential associated to some of these charges. If the chemical potential is sufficiently high one can 
show that the gaps (i.e. the energy at zero momentum) of these vectors become light [OJ [l2j and eventually zero 
signaling an instability. If one applies our results to 2 color Quantum Chromo Dynamics (QCD) at non zero baryon 
chemical potential one predicts that the vectors made out of two quarks (in the 2 color theory the baryonic degrees 
of freedom are bosons) condense. Recently lattice studies for 2 color at high baryonic potential jf3| seem to support 
our predictions. This is the relativistic vectorial Bose-Einstein condensation phenomenon. A decrease in the gap of 
vectors is also suggested at high baryon chemical potential for two colors in [ p~i| . 

Interestingly non-relativistic vectorial Bose-Einstein condensation recently has attracted much attention in con¬ 
densed matter physics since it has been observed experimentally in alkali atom gases |h|. We also note that in this 
framework a rich phenomenology related to the classical solutions of the theory, such as vortices (l(| , is expected to 
occur. 

Here we extend the analysis presented in 00 by showing the existence of new phases while providing a detailed 
investigation of the relativistic massive vector condensation phenomenon driven by a non zero chemical potential. 
We develop an efficient formalism which allows us to analytically study the dispersion relations of the vectors in the 
broken phase of the theory and set the stage for possible higher order computations. Our Lagrangian approach has 
to be understood as a Landau theory describing the vector condensation phenomenon. 

Having in mind some specific physical applications we assume our massive relativistic vectors to be in the adjoint 
representation of the non abelian SU(2 ) global symmetry group. We then turn on a chemical potential in a specific 
SU( 2) direction which breaks SU(2) explicitly to a 17(1) symmetry. We demonstrate that we can have three inde¬ 
pendent phases according to the value assumed by one of the two distinct (non derivative) vector self interactions. 
The polar phase, the apolar phase and the enhanced symmetry one. The polar phase has been introduced in j[T^1 
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along with the enhanced symmetry case. The polar and apolar phases are the relativistic generalization of the phases 
encountered in the condense matter |l^] framework. The appearance of a given phase is related to the specific value 
of one of the vector self interaction coefficients of the Lagrangian at zero chemical potential. In this paper we study 
consistently the gapless excitations and the dispersion relations in all of the phases while providing a formalism helpful 
when trying to go beyond the tree level approximation. 

In the polar phase the vector condensate breaks the rotational symmetry down to a simple rotation around the axis 
of condensation. The internal U{ 1) symmetry breaks completely as well. This phase is characterized by three gapless 
excitations with linear in momentum dispersion relations. By studying in detail the dispersion relations we show that 
two of the gapless states have isotropic dispersion relations while the third excitation has different velocities in the 
direction parallel and orthogonal to the vacuum expectation value. The detailed analysis of the dispersion relations is 
exact when we supplement the global symmetry SU (2) by an extra one which allows only an even number of vectors 
in any vertex of the theory, however the number and type of goldstone bosons is independent of the extra discrete 
symmetry. 

If we are in the apolar phase the condensation of the vector is such that the unbroken generator is a linear 
combination of one of the rotational generators and the abelian 17(1) generator. Only two gapless excitations emerge 
this time. One with linear and the other with quadratic dispersion relations. However the vacuum still breaks three 
generators. This is in agreement with the Nielsen and Chadha 0 counting scheme according to which in absence of 
Lorentz invariance (and under a number of assumptions) each gapless excitation with quadratic dispersion relations 
has to be counted as two goldstone bosons with linear dispersion relations. More specifically if nj denotes the number 
of gapless excitations of type I with linear dispersion relations (i.e. E oc p) and nj/ the ones with quadratic dispersion 
relations (i.e. E oc p 2 ) the Goldstone theorem generalizes to nj + 2njj > # broken generators. Clearly in the absence 
of type II excitations we recover the usual counting. We analyze the spectrum in this case and find that the type 
II goldstone boson has isotropic dispersion relations while the type I has not. In all of the phases the non goldstone 
boson excitations are investigated in detail as well. 

If one of the self-vector couplings is set to zero we observe three gapless excitations: 2 type II goldstone bosons and 
one type I. In this case the potential has an enhanced global symmetry, i.e. it has an SO( 6) which the vacuum breaks 
to an SO( 5). In absence of Lorentz breaking, we would have 5 ordinary goldstone bosons. However the presence of 
the chemical potential in the derivative terms of the Lagrangian prevents the emergence of 2 extra goldstone bosons 
]l8l while turning 2 goldstones into type II ]l2| |. 

We also discover that when the vector self interaction couplings are tuned to be the ones predicted using the Yang- 
Mills relation (while keeping always a non zero vector mass) the theory at high chemical potential does not predict a 
stable solution. The existence of an inhomogeneous phase is not excluded. However a more natural solution of the 
instability is that the chemical potential, when increasing the associated charge density, is at the most as large as the 
mass of the vectors. 

If we now apply our results to describe non elementary relativistic massive vectors such as the p field of QCD at high 
chemical potential (isospin for example) it will help elucidating and constructing the effective Lagrangians for vector 
fields at zero chemical potential. Interestingly, indeed, the particular choice of the vector coupling respecting Yang- 
Mills relations is the one often assumed in literature when introducing the quark composite field p at the effective 
Lagrangian level jl9j ^o|. For example the gauge choice of the couplings emerges very naturally when the vector 
mesons are introduced as gauge fields of an hidden local gauge symmetry js) . By studying the vector condensation 
phenomenon (for example the ordinary p at high isospin chemical potential) on the lattice it is possible to shed light on 
the correct way of constructing a theory for composite vector fields or massive vector fields in general. These effective 
theories for composite fields are relevant also for the physics beyond the standard model of particle interactions such 
as the ones relative to a strongly interacting electroweak sector (see HI-IH). 

In section |l| we present the Lagrangian and study the vacuum. We show that a number of phases can be present. 
In sec. 0 we solve for the dispersion relations. We then conclude while reviewing the physical applications in sec. |Tv| . 
In the appendices we summarize first all of the terms of the Lagrangian in the cartesian and cylindrical coordinates. 
We finally show the computational details of the dispersion relations for the apolar case. 


II. VACUUM STRUCTURE AND DIFFERENT PHASES 

There are different ways to describe vector fields at the effective Lagrangian level. For example one can use the 
hidden local gauge symmetry of Ref. Q, or the antisymmetric tensor field of Ref. ^3|. Or one can introduce the 
massive vector fields following the method outlined in 00- For most of the known physical applications these 
approaches provide identical results (at the tree level). It is worth noticing that in some of these methods certain 
coefficients of the vector Lagrangian are related by enforcing the Yang-Mills relations. 

We choose to consider the following general effective Lagrangian for a relativistic massive vector field in the adjoint 
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of SU ( 2 ) in 3 + 1 dimensions and up to four vector fields, two derivatives and containing only intrinsic positive parity 
terms |26j : 


C = --F^F a ^ + 


nn 


+ 6 e abc 


dpA c 


y A^A:-j(A;A ^) 2 + j(A;A^y 


(1) 


with = dpA'p — d l ,A“, a = 1, 2, 3 and metric convention 77 ^ = diag(+, —, —, —). Here, <5 is a real dimensionless 
coefficient, m 2 is the tree level mass term and A and A' are positive dimensionless coefficients with A > A' when 
A' > 0 or A > 0 when A' > 0 to insure positivity of the potential. The Lagrangian describes a self interacting SU (2) 
Yang-Mills theory in the limit m 2 = 0, A = A' > 0 and <5 = — \f\. 

It is relevant to notice that in the limit <5 = 0 the theory gains a new symmetry according to which we have always 
a total number of even vectors in any process. This symmetry guarantees that if the <5 term is absent from the start 
it will not be generated dynamically. In this paper we will mainly investigate the theory in this case since it will 
simplify our computations. However we will comment on the effects of such a term in a final paragraph. 

r 3 

The effect of a nonzero chemical potential associated to a given conserved charge - (say T 3 = —) - can be readily 
included [pT| by modifying the derivatives acting on the vector fields: 


d u A p -> d v A p — i [B v , A p ] 


( 2 ) 


with B v = fiSvoT 3 = V v T 3 where V = (n ,0). In appendix it is summarized the effective Lagrangian, in this basis, 
after the introduction of the chemical potential. The introduction of the chemical potential breaks explicitly the 
Lorentz transformation leaving invariant the rotational symmetry. Also the SU (2) internal symmetry breaks to a 
1/(1) symmetry. If the 6 term is absent we have an extra unbroken Z 2 symmetry which acts according to A 3 —> — A 3 . 
These symmetries suggest introducing the following cylindric coordinates: 

</V = ^(K + iAl ) , = -^=(A 3 - 5A 2 ) , % = A 3 , (3) 

on which the covariant derivative acts as follows: 


Dp,<t> v = (d+ iV) fi <j> I/ , Dpip v = dpil) v , V v = {n, 0 ) . (4) 

The quadratic, cubic and quartic terms - in the vector fields - in the cylindrical coordinates are summarized in the 
appendices. 


A. The Non Derivative Terms 


We first study the non derivative 

V = —~y [20*V° + ip •V'] + 

■ il>) - 


terms which we collect in the following potential type term: 


9 \ _ \/ y 

( m 2 - / i 2 ) [</>* V ° -(/)*■ ( f >\ +--- { 4 >* ■ ( j )) 2 - —1(0 • <£)| 2 + 

■ 4>)\ 2 + 2[i5 %lP (4>* ■ 0) + /j,S (j)*° (cj) ■ tp) + c.c. 


A - A '// M2 

—4— (V’-VO 


( 5 ) 


We can read off the symmetries of the theory from the potential. When <5 = 0, for example, we gain the discrete 
symmetry ij} —> —ip. To explore the vacuum structure of the theory we consider the following variational ansatz: 


r = 0 , 


= <T 



( 6 ) 


Substituting the ansatz in the potential expression we have: 

V = 2 a 4 [(2A - A') - A' cos 2 a] +2 (m 2 - ^ 2 ) cr 2 . (7) 

The potential is positive for any value of a when A > A' if X' > 0 or A > 0 if A' < 0. Due to our ansatz the ground 
state is independent of <5. The unbroken phase occurs when fj, < m and the minimum is at cr = 0. A possible broken 
phase is achieved when /i > m since in this case the quadratic term in a is negative. According to the value of X' we 
distinguish three distinct phases: 
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B. The polar phase: A' > 0 


In this phase the minimum is for 

/ 

= o' 

V 

This phase has been partially analyzed in [^2| . 

SO(2). If we define with Ri the three generators of SO (3) (acting only on the spatial indices) as follows: 

/ 0 1 o\ / 0 0 1\ / 0 0 0 \ 

Ri = -1 o 0 , R 2 = 0 0 0 , R 3 = 0 0 1 , (9) 

\ o o o/ y-i o o/ yo-ioy 

while the 17(1) generator acts as a phase on <f> then the unbroken generator determined imposing invariance of the 

vacuum 


0 


11 1 I /" _ 1T| " 

| , with = , 


4 A-A' 

We have the following pattern of symmetry breaking SO( 3) x 17(1) 


( 8 ) 


^untokenW = 0 , 


( 10 ) 


is the linear combination T^° b ” ken = R 2 — f? 3 - We have three broken generators. We can show (see section on the 
dispersion relations for details) that in this case we have 3 gapless excitations with linear dispersion relations. All of 
the physical states (with and without a gap) are either vectors (2-component) or scalars with respect to the unbroken 
50(2) group. The dispersion relations for the 3 gapless states Q computed explicitly in the following sections are 
reported here for the reader’s convenience. 


E l- 


oc A (/u — to 2 ) p , 


V = 2 Physical States 


Eg - oc 


(/i 2 — m 2 ) + n 2 s _p|^ , S = 1 Physical State 


« (A-A') 


, fi 2 A + A' (/i 2 — 2m 2 ) 


(p 2 A - m 2 A ') 2 


( 11 ) 

( 12 ) 

(13) 


where P||(±) refers to the momentum parallel (orthogonal) to the vacuum. Here we present the leading terms in a 
momentum expansion of the gapless dispersion relations. The vector components orthogonal to the vacuum direction 
(2 states indicated with V) propagate isotropically while the component in the direction of the vacuum (1 state 
indicated with S) does not. As a consistency check one sees that at /i = in the dispersion relations are all isotropic. 
At [i = m the dispersion relations are no longer linear in the momentum. This is related to the fact that the specific 
part of the potential term has a partial conformal symmetry discussed first in JT^] . Some states in the theory are 
curvatureless but the chemical potential term present in the derivative term prevents these states to be gapless. There 
is a transfer of the conformal symmetry information from the potential term to the vanishing of the velocity of the 
gapless excitations related to the would be gapless states. This conversion is due to the linear time-derivative term 
induced by the presence of the chemical potential term II @- 


C. Enhanced symmetry and type II Goldstone bosons: X' = 0 

Here the potential has an enhanced 50(6) in contrast to the 51/(2) x 17(1) for A' ^ 0 global symmetry which breaks 
to an 50(5) with 5 broken generators. Expanding the potential around the vacuum we find 5 null curvatures [ |l2| . 
However we have only three gapless states obtained diagonalizing the quadratic kinetic term and the potential term 
(see the dispersion relation section). Two states (a vector of 50(2)) become type II goldstone bosons (see eq. @) 
while the scalar state remains type I. This latter state is the goldstone boson related to the spontaneously broken 
U( 1) symmetry. According to the Nielsen-Chada theorem the type II states are counted twice with respect to the 
number of broken generators while the linear just once recovering the number of generators broken by the vacuum. 
More specifically one can prove [|l2j that the velocity of the states labelled by V in eq. ( 0 ) is proportional to the 
curvatures (evaluate on the minimum) of the would be goldstone bosons which is zero in the X' = 0 limit. Again we 
have an efficient mechanism for communicating the information of the extra broken symmetries from the curvatures 
to the velocities of the already gapless excitations. 
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D. The apolar phase: X 1 < 0 


We now extend the vacuum analysis in 


0 


to the apolar phase in which the potential is minimized for: 


(r) 



with a 2 


1 /x 2 — m 2 

2 2A - A' 


(14) 


In this phase we have again 3 broken generators. However the unbroken generator is the following combination of the 
f/(l) and SO( 3) generator: 


^unbroken = ^1 “ * * ■ (15) 

More explicitly the action on the vev is (Ri — 1 1 )(cf>) = 0. In this phase as we shall demonstrate in the next sections 
we have only two gapless states. One of the two states is a type I goldstone boson while the other is type II. The two 
goldstone bosons are one in the 2 and the other in the x — y plane. Interestingly in this phase, due to the intrinsic 
complex nature of the vev, we have spontaneous CP breaking. 

We summarize in Fig. jll the phase structure in terms of the number of goldstone bosons and their type according 
to the values assumed by A'. 


1 type I 
1 type II 


t 


3 type I 


A' 


1 type I 

2 type II 


FIG. 1: We show the number and type of goldstone bosons in the three distinct phases associated to the value assumed by 
the coupling X'. In the polar phase, positive X ', we have 3 type I goldstone bosons. In the apolar phase, negative A', we have 
one type I and one type II goldstone boson while in the enhanced symmetry case A' = 0 we have one type I and two type II 
excitations. 


E. The case A = Ah the gauge theory limit 


Here the potential is: 

V = 2 ex 4 A sin 2 a + 2 ( m 2 — p 2 ) a 2 . (16) 

This potential has two extrema when /i > to, one for a = 0 and a = 0 which is an unstable point and the other 
for a = ±7t/ 2 and a 2 = ^ jjr corresponding to a saddle point (see the potential in Fig. |^). At first the fact that 
we have no stable solutions seems unreasonable since we know that in literature we often encounter condensation 
of intermediate vector mesons such as the W boson. However (except for extending the theory in higher space 
dimensions) in these cases one often introduces an external source. For example one adds to the theory a strong 
magnetic field (say in the direction z) which couples to the electromagnetically charged intermediate vector bosons 
W + and W~). In this case the potential is (see ||): 

V = 2 er 4 A sin 2 a + 2 (to 2 — e H sin a) a 2 , (17) 

where e is the electromagnetic charge and H is the external electromagnetic source field. This potential has a true 
minimum for a = 7r/2 and a = e 2 a™ whenever the external magnetic field satisfies the relation eH > to. We 
learn that the relativistic vector theory is unstable at large chemical potential whenever the non derivative vector 
self interactions are tuned to be identical. This is precisely the limit often used in literature when writing effective 
Lagrangians that in QCD describe the p vector field. In principle we can still imagine to stabilize the potential in the 
gauge limit by adding some higher order operators which seems unnatural. A more natural solution to this instability 
is that the chemical potential actually does not rise above the mass of the vectors even if we increase the relative 
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charge density. This phenomenon is similar to what happens in the case of an ideal bose gas at high chemical potential 

& 

Interestingly by studying the vector condensation phenomenon for strongly interacting theories on the lattice at 
high isospin chemical potential we can determine the best way of describing the ordinary vector self-interactions at 
zero chemical potential. 


III. DISPERSION RELATIONS 


To determine the dispersion relations we concentrate on the quadratic terms of the theory. For n < m vectors do 
not condense and the only terms we need to consider are the ones in eq. (B2) which we report here for the reader’s 
convenience: 


^Quadratic = \^ \Sfiv [Q 2 + TO. 2 ] - d^d„] Ip" + j ^ (g^[D 2 + TO 2 ] - D^D V ) 4? + C.C.j . 


(18) 


These are the only terms we need for the SU( 2) theory with or without the S term. We note that in the covariant 
derivative acting on <fi it is hidden the negative /i 2 square term appeared already in the potential term in the previous 
section. The ip held is a standard massive vector held with 3 independent degrees of freedom since it satishes the 
constraint m 2 d ll ip IJ ' = 0. The associated dispersion relations for the 3 physical components of the ip fields are: 

E^j = y/p 2 + TO 2 . (19) 


For the field (p we have the following equation of motion 

(<fo„ [D 2 + m 2 ] - Dfj,D v ) cf> v = 0 , (20) 

which by multiplying on the left by D M leads to the constraint m 2 D l/ <j> 1 ' = 0 reducing the number of physical degrees 
of freedom for to three. Assuming the previous constraint the equation of motion is clearly [D 2 + m 2 ]<p v = 0, 
leading to the following dispersion relations jf2| : 

E = ±fj, + y/p 2 + m 2 . (21) 

Each of the (±) state corresponds to a positive and negative charge under 1/(1) and constitutes an SO(3) vector with 
3 independent degrees of freedom. Having reviewed the case g < m we now analyze the dispersion relations of the 
system when the vector condenses (i.e. /i > to). The physical constraints are now more involved. Expanding our 
fields around the new vacuum of the theory: 


<t>li = (0m) + 0m > 0M = 0m > ( 22 ) 

new quadratic terms emerge depending on (0 M ). Since the specific properties of the vacuum are very different for 
A' > 0 or A' < 0 we consider these two cases separately. 
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A. The polar phase dispersion relations 

1. Dispersion relations for the if field 


In the 5 = 0 limit the quadratic terms of if and <f decouple and we start analyzing the dispersion relations for the 
field if. The most general Lagrangian term quadratic in the fluctuation field if is: 

^ [g^i® 2 + m 2 } - 9fj,d„] if v + ^if^n^,n v if 1 ' , (23) 

with = {(ffi) and 

to 2 = m 2 — 2 Xn^ri/j, , k = 2 A' . (24) 

Here we have already used the fact that is real. Multiplying the associated equation of motions on the left by 
the free field constraint is: 


m 2 d v if v = —k n ■ d n v if 1 ' 


(25) 


A convenient way of dealing with this constraint is to split our field in a component parallel and one orthogonal to 
the vacuum expectation value via [E7[: 


/ -r> if i n n n f iv _ / i n n 

^Vu H- — 


■;^S , 


with 


V = - 

’ uf y 




n^n 


V-r 


tfvu = 'P^if v , if s = 


n u if v 

\/—n 2 


(26) 


(27) 


Clearly ify is transverse to the vev while ifs is parallel. In terms of these new fields and using the constraint in 
eq. (Eq) the quadratic Lagrangian for if as function of a generic but real vev for cf is: 


1 


-ifv^ ( d 2 + m 2 ) if Vll - -if s 


( d 2 + to 2 ) + n 2 k - n 2 — (n ■ d ) 2 


ifs ■ 


For the 2-conrponents transverse field ify we derive the following dispersion relations: 


e 2 v = p 2 + A 2 v 


= ™ 2 


while for the (one component) longitudinal field ^$ we have: 

td 2 2,2 2 1 a 2 a 2 -2,2 

E 1>s =P±+ v l>s\\P II + ’ A bs = TO + n k 

Substituting the explicit expression for the vacuum we deduce: 


l Vs II 




m 


X 


A lv =P 2 + T—T7 (m 2 - m2 ) > A ls = P 2 > 


A - A' vr ~ ' ’ _vs " ’ M - (A-A') 2 m 2 

P\\(±) are the component of the momentum parallel (perpendicular) to the vector condensate. 


= 1 - 


A' 


/Z — TO 


(28) 

(29) 

(30) 

(31) 


2. Dispersion relations for the <f field 


The situation is more involved for the fluctuations of the complex <f> vector field. For these fields we write the 
general quadratic Lagrangian in a 2 component formalism as follows: 

TjX^ [_9nf (A 2 + M) - A m A„] x u + \x ]li ■ (32) 

with 


X = 



A. 


Dn 0 \ 

o d;J 




a 

p 


(33) 












The equation of motion for \ is: 


[fiV (A 2 + M) - A^ A„] x v + ICn l _ l n„x 1 ' = 0 • 

Multiplying on the left by we obtain the following physical condition: 


MA„x v = - (n • A) /C (n • x) • 


We again split Xu as follows: 


*v IJ: = 'PfivX 1 ' 


and using eq. (pq) the quadratic Lagrangian becomes: 


= 


n„X 

y/—n 2 


i$£ f (A 2 + M)$ V(U -(A 2 + M) + n 2 /C - n 2 K ] M X M l K{n-df 




(34) 

(35) 

(36) 

(37) 


In deriving the last equation we used the fact that the zeroth component of vanishes. When the covariant derivative 
acts on $ is always D while when acts on $* is always D*. We have an expression, formally, similar to the one obtained 
for the neutral field i/j. 

In terms of the coefficients of the effective Lagrangian we have (after substituting the expression for the vev): 


A 2 = 

K = - 


-£ 2 + p 2 + 2/t£^f( 2 

° ^ 

0 

-E 2 +p 2 -2/iE- [i 2 ) 

'2A-3A' 2A — A' \ 

2 1 fj, 2 — m 2 

v 2A - A' 2A — 3A' ) ' 

’ U 2 A-A' ‘ 

case holds the relation: 

—2 2 

nn — (i = 


(m 2 -( 2A-A') n 2 
M - 1 X'n 2 


f~,2 


A ’n 

to 2 — (2A — A') n 2 


(38) 


(39) 

we have one zero gap vector field (with 2 independent physical degrees of freedom) with respect to 50(2) and one 
non zero gap vector with 


A 2 + = 2 (to 2 + /r 2 ) = 4/r 2 + A' 


fji 2 — TO 2 

A - A' 


The full dispersion relations are (without enforcing yet the condition (|39|)): 


E 2 ^ T = m 2 + p 2 - p 2 =F V 4/r 2 p 2 + r 2 + 4?n 2 
Expanding in momenta and enforcing eq. 

A' n 2 — m 2 


Ei- = 


A - A' A 2 


P 2 + 16 Xf^P 4 + °^ 6 ^o") ’ 


A 6 


$ 


+ 




+ 


E : 


;+ = + 1 


4 ii 2 

AE 


P 2 -16^P 4 + 0(m 6 ^), 


*7 


All 


(40) 

(41) 

(42) 

(43) 


*7 


Hence the propagation of the vector fields orthogonal to the vev is isotropic. 

We are now left to investigate >I>s. Using eq. ([f7|) we get the following dispersion relations for the gapless mode 


9 2 

/i — m 


E 2 = 

45 3 [i 2 — rri 


— 5 [Pl + vl-pU+--- , 


v l s - = M 2 (A A' 


M ^ 2 A + A' (n 2 - 2m 2 ) 

(/x 2 A — m 2 A') 2 


( 44 ) 


The massive mode dispersion relation is: 


El s+ = 2(3 M 2 -m 2 ) + 


5/i 2 — m 2 


3/z 2 — m 2 


2 P±+^ s +P|| +•" > 
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with vi + a lengthy but known expression of the Lagrangian coefficients. This completes the analytical study of the 

®S|| 

dispersion relations for the case of the polar phase. We learn that the vector states orthogonal to the condensate have 
isotropic dispersion relations while the ones relative to the vector component in the direction of the condensate are 
not isotropic. We also find that for X' = 0, as anticipated in the previous section, 2 gapless excitations have quadratic 
dispersion relations and hence become type II goldstone bosons. This is related to the enhancement of the global 
symmetry in the potential term. In Fig. [3| we plot the gaps as function of the chemical potential for this phase in the 
left panel. Before condensation each solid line corresponds to three physical states, after condensation three gapless 
modes emerge and each dashed line corresponds to two states while each solid one to a single state. 




fi/m jt/m 


FIG. 3: We present the Gaps ( E(p = 0)) for the polar phase (left panel) and the apolar (right panel) of the theory as function 
of the chemical potential in units of m. Left Panel: Before condensation each line describes three degenerate massive states. 
After condensation fj, > m the dashed lines represent two states while the solid lines represent one state each except for the 
three gapless states. We used the following values for the plot: A = 1, A' = 0.33. Right Panel: Before condensation each line 
describes three degenerate massive states. After condensation fj, > m the dashed line represents two states (ipv) while the solid 
lines one physical state each except for the two gapless states. We used the following values for the plot: A = 1, A' = —0.3. 


B. The Apolar Phase Dispersion Relations and Gaps 

We also solved the physical constraint for the apolar case. In this case the analytical analysis of the physical con¬ 
straints and of the dispersion relations is complicated by the fact that the vacuum is complex. Since the computations 
are instructive but technical we provided them in the last appendix and summarize here the results. In this phase we 
have three broken generators but only two gapless modes: a type I and a type II goldstone boson. Since the unbroken 
generator is a linear combination of a rotation and the generator for the internal 17(1) symmetry the gaps and the 
dispersion relations lose the straightforward and nice classification in doublets (i.e. vectors) and scalars with respect 
to a standard rotation. The gaps in this phase, and for a specific choice of the couplings, are displayed in the right 
panel of Fig. [| Before condensation each line describes three degenerate massive states. After condensation /i > m 
the dashed line (ipv) represents two states while the solid lines one physical state each except for the two gapless 
states. We used the following values for the plot: A = 1, A' = —0.3. The dispersion relations too, in the <5 = 0 limit, 
display a more complex structure which does not alter the gapless excitation structure and goldstone counting. 


C. Breaking the Z 2 symmetry. 

We now comment on what happens if we allow Z 2 symmetry breaking terms such as the S term. Before including 
the chemical potential in the direction T 3 the absence of the <5 term prevents an odd number of vectors to be present in 
any vertex of the theory. Since this term involves 3 fields it will not affect the dispersion relations before condensation. 
After condensation has taken place the vacuum structure (due to our ansatz) is also unaffected by this term. Since 
the goldstone states are the fluctuations around the vacuum in the direction of some of the continuously broken 
symmetries their general properties are also expected not to be disrupted. The S term will, however, change some 
of the details of the dispersion relations. The possible quadratic terms in the fields emerging after condensation will 
always mix the %l> state with a (j> one and their effect will be investigated elsewhere. 
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IV. PHYSICAL APPLICATIONS AND CONCLUSIONS 

We investigated the phase structure of the relativistic massive vector condensation phenomenon due to a non zero 
chemical potential associated to some of the global conserved charges of the theory. The possible phase structure is 
very rich. Indeed according to the value assumed by X 1 we have three independent phases. The polar phase with X 1 
positive is characterized by a real vacuum expectation value and 3 goldstone bosons of type I. The apolar phase for 
A' negative has a complex vector vacuum expectation value spontaneously breaking CP. In this phase we have one 
goldstone boson of type I and one of type II while still breaking 3 continuous symmetries. The third phase has an 
enhanced potential type symmetry and 3 goldstone bosons one of type I and two of type II. 

We also discovered that if we force the self interaction couplings A and A' to be identical, as predicted in a Yang- 
Mills massive theory, our ansatz for the vacuum does not lead to a stable minimum when increasing the chemical 
potential above the mass of the vectors. A possible resolution of such an instability is that the chemical potential can 
be at most as large as the vector mass. This case it very similar to the Bose-Einstein condensation phenomenon for 
an ideal bose gas at high chemical potential. Interestingly the gauge coupling limit is often adopted in literature to 
economically describe, for example, the QCD composite vector field p. We hence suggest that lattice studies at high 
isospin chemical potential in the vector channel for QCD might be able to, indirectly, shed light on this sector of the 
theory at zero chemical potential. More generally the hope is that these studies might help understanding how to 
construct consistent theories of interacting massive higher spin fields not necessarily related to a gauge principle. 

We also developed a formalism which enabled us to investigate the vacuum structure and dispersion relations in 
the spontaneously broken phase of the theory. Our results are helpful when trying to go beyond the classical and tree 
approximation. Our present studies are readily applicable to a number of physical phenomena of topical interest. For 
example in the framework of 2 color QCD at high baryon chemical potential vector condensation has been predicted 
in 00- Recent lattice studies |l3| seem to support it. The present analysis while reinforcing the scenario of vector 
condensation shows that we can have many different types of condensations with very distinct signatures. Some 
details of other possible physical applications are presented in Q. The present analysis can be straightforwardly 
extended to a general number of space dimensions ||l2| which may be useful for more exotic scenarios related to the 
phenomenon of vector condensation Q. 
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APPENDIX A: LAGRANGIAN IN CARTESIAN COMPONENTS 


It is helpful to know also the different terms of the theory in cartesian components. We start with the vector 
Lagrangian presented and studied in the main test: 


1 777 ^ 

C = --F^F^F — A“A ap + 6e c 


,bc 


dpA a „A»A: - ± (a;a^) 2 + j {a;a^Y 


(Al) 


with = d^A® — a = 1, 2, 3 and metric convention rj^ = diag(+, —, —, —). Here, 6 is a real dimensionless 

coefficient, m 2 is the tree level mass term and A and X 1 are positive dimensionless coefficients with A > X'. 

After including a nonzero chemical potential associated to a given conserved charge - related to the generator (say 
B) - in the the following way: 

d v A p -> d v A p — i [B„ , A p ] , (A2) 

with = pl5 v qB = V„B where V = (p, ,0). The vector kinetic term modifies according to: 


Tr \Fp V F pv ] - Tr[F p „F pI '}-4iTr[Fp V [B p ,A I '])-2Tr[[B p ,A v }[B p ,A 1 '}~[B p ,A v ][B l ',A p }} . 
The terms induced by F^F^, after integration by parts, yields [II]] 


(A3) 


= \a% l S ab [g pv U - d p d v ] - 4i 7ab 


g pv V ■ d - 


V p d v + V v d p 


+ 2Xab [V ■ Vg pv - V»V V ] \ A 


kinetic 


(A4) 
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For B = T 3 we have 


lab = Tr [T a [B, T b ] ] , Xab = Tr [[B, T a } [B, T b ] . 


X 1 

lab = ^2 £ah3 ’ Xll = X22 = — g ) X33 = 0 ■ 


The chemical potential induces a “magnetic-type” mass term for the vectors at tree-level. 
The trilinear term with a single derivative is: 

C s = 5e abc d^A av A^A" + 5 (V ■ A b ) A a " A c „ [<5 36 <5 oc - <5 3c <5 ab ] . 


APPENDIX B: CYLINDRICAL COORDINATES 


Here we summarize all of the terms in the Lagrangian using the cylindrical coordinates: 

K = ±(Al + iAl) , <p*» = ±(Al - iAl) , % = Al , 

The quadratic, cubic and quartic terms - in the vector fields - now reads: 

^quadratic = \^ [SV[<9 2 + m 2 ] - <9^] V’ 1 ' + (p* 11 {9 iiu[D 2 + m 2 ] - D^D V ) <p" + c.c. 


Ccubic = -i 0* M • <p") - i Sip" (d^cp*^) - i Sip* 1 {d li <j>l<p v ) + c.c. 
—25 V ■ ip ((p* ■ cp) — 5 [V ■ <p* (<p ■ ip) + c.c.] . 


Creators = ~(A;A^) 2 + ^(A;A^) 2 = (0* ■ <Pf + ^ | (0 ■ *) | 2 + ^ & ■ V-) 2 


- A (</>*•</>) ( V > • VO + A ' | 


APPENDIX C: DISPERSION RELATIONS FOR THE APOLAR PHASE 


In this case the vev is proportional to the vector: 


/° 

J_ 1 

V2 ' 1 


with n 2 = n* 2 = 0 , and n ■ n* = —1 


The quadratic Lagrangian term for the held ip takes the form: 


\ V [9^[d 2 + m 2 ] - d^dv] ip" + (n M n* + n*n„) ip" , 


and the constraint equation reads: 


m 2 d ■ ip = —k [n ■ d (n* ■ ip) + n* ■ d (n ■ ip)] . 


We define the following projectors 


Pfill 9flU PflV ) P fill 


n^nl + n„n* 
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Now has two non zero components (i.e. az = v = x and /i = v = y) and hence projects out two real scalars while 
has non zero the temporal and the zed component projecting out a vector which lives in 1 + 1 dimensions (i.e. 
again a scalar field). It is convenient to split z/> as follows: 


VV = ^Sn + 


Tin 


-ipD 




n • n* n • n* 

Using these fields plus the constraint the Lagrangian becomes: 


D > 


ip D = n* ■ tp , 


£ = (<9 2 + to 2 ) - ip* D [d 2 + TO 2 + ^-j(n • <9)(n* • <9)^ ip D 

“2^4 [f/ ) D(n-a) 2 fe+ft(n*'9) 2 ^] , 


with 


A 2 _ ~ 2 —2 

A 0d = m =m ~ K > 


A 2 s = m 2 . 


Diagonalizing the D and S' sector independently we deduce the following dispersion relations: 


= a I.+p 2 


with 


Ud 2 — A V>d + V D2 (Px + Py) + Pz ) — TO 2 + P 


= 1 - ^4 • 


Substituting for the vev: 


= TO = At 


A' 


2A-A' 


(a t 2 — m 2 ) , k 


A' 


2A-A' 


(p 2 — m 2 ) , A 2 D = az 2 . 


(C5) 


(C 6 ) 

(C7) 

(C 8 ) 

(C9) 

(CIO) 


We see that in the z/> sector of the theory when S is zero two states have isotropic dispersion relations but are not 
degenerate. The two degenerate states have different momentum dependence. 


1. The <f> dispersion relations 

The general quadratic Lagrangian in a 2 component formalism and in the presence of the complex vev for <j> reads: 

[AW (A 2 + M) - A„A v ] x" + l -x^^ vX v ■ (Cll) 

with 

AZ 2 0 


M = 


0 a* 


, = (to 2 - /. L ) 


n n n v — 'yn^riv 

~ 7 WX 


3 k >: 
n ^ n V 


2A' 


7 = 


2A-A' 


(C12) 


Note that the matrix Lorentz structure is solely determined by the vector vev. The quadratic equation of motion 
leads to the following physical constraint: 


MKx v = -A ^ vX v 


We split x as follows: 


Xfj. — + 


n u ~ n u 

- -&d -\ --d>D , 


= n* • x , <E> D = n • x • 


(C13) 


(C14) 


Now $ 5 , d>£> and d>£> represent three independent (2-components) physical states, (i.e. using these fields plus the 
constraint the quadratic term Lagrangian can be compactly written as: 


Ugt (A 2 + M) - \r D t 


(A 2 + M)xl + Uq 


Td , 


(C15) 
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with 


( \ tt — ( n *lOddK,n n*lOd dJCn* \ „ _ / n*ICn n*lCn* \ 

V / ’ U ~ V ntfddKn nK)ddK,n* ) ’ Q ~ ^ nICn nICn* J 


(C16) 


where Tp is a 4 column vector, U and Q are four-dimensional matrices. nlOd = n ti lO flu d 1 ' is a two dimensional 
matrix. 

The dispersion relations in the S sector are now straightforward and we deduce the following two states: 


Egb-ii = £ 2 + <V) , (C17) 

Ewbgb = V + |^ + 0(p A ) . (CIS) 

The first state is a type II goldstone boson while the second state is the would be goldstone boson e.g. the one which 
would have been massless if we had no breaking of the Lorentz symmetry. For the D sector the diagonalization can 
be performed analytically when setting the momentum in the x and y direction to zero and we get: 

E 2 G b-i = v z 2 p 2 z + 0(p 4 ) , (C19) 

E% a = 2(3p 2 - m 2 ) + (2 + vl)p 2 z + 0(p 4 ) , (C20) 

E\ b = A| b +u| b pI + 0(p 4 ), (C21) 

E% c = + v% c p 2 + 0{p 4 ) . (C22) 


This sector of the theory contains a goldstone boson of type I and 3 massive states. 
In the following we summarize the nine physical gaps related to this phase: 


A 2 

^GLess 

= 0 , 

2 States 

^ <f>WBGB 

= 4 p 2 

1 State 

a 5„ 

= 2 (3 p 2 — m 2 ) 

1 State 

a 5„ 

~ (2A-A /) 2 (^-1 + V^2) 

1 State 

O 

(N-O- 

<1 

— (2A-A') 2 (^ Vh2) 

1 State 


(C23) 


A !v = M 2 

2 States 

A L = P 2 + ^V(P 2 -™ 2 ) 

1 State 


(C24) 


with: 


hi = m 2 { 2A - A')A' + 4A 2 /i 2 - 6\X'p 2 + 2 A'V, h2 = p 2 ( 2 A - A ') 3 (2m 2 A' + (2A - 3A ')p 2 ) . (C25) 
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